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In  the  early  1950's  the  Corps  of  Engineers'  Jacksonville  District 
initiated  a series  of  laboratory  tests  to  investigate  the  overtopping  of 
proposed  levee  sections  for  Lake  Okeechobee,  Florida.  For  economic 
reasons,  the  alternative  to  build  levees  with  crest  elevations  that  were 
at  times  below  the  limit  of  wave  runup  was  investigated  and  the  quantities 
of  water  carried  over  the  structures  for  various  freeboard  allowances, 
structure  slopes  and  wave  conditions  determined.  The  initial  tests  were 
conducted  at  the  Waterways  Experiment  Station  (WES)  in  Vicksburg, 
Mississippi  for  the  Jacksonville  District  at  what  was  taken  to  be  a 
1 to  30  model  scale.  Model  wave  heights  varied  from  I4.05  cm  to  12,2  cm 
(0,133  to  O.UO  ft).  In  order  to  expand  the  range  of  test  conditions 
investigated,  the  Beach  Erosion  Board,  currently  the  Coastal  Engineering 
Research  Center  (CERC),  commissioned  an  expanded  series  of  tests  that 
considered  the  overtopping  of  riprap  faced,  curved  and  stepped  seawalls 
as  well  as  the  overtopping  of  "smooth"  slopes.  These  tests,  also  con- 
ducted at  WES,  were  considered  to  be  at  a 1 to  1?  scale  with  model  wave 
heights  ranging  from  5*36  cm  to  21.5  cm  (0.176  to  0.706  ft).  A number 
of  tests  were  subsequently  conducted  in  CERC 'a  large  wave  tank  to 
determine  the  influence  scale  effects  might  have  on  overtopping.  These 
tests  are  referred  to  as  1 to  2 1/2  scale  tests.  The  model  wave  heights 
investigated  ranged  from  i<8.0  cm  to  1^0.2  cm.  (I.60  to  h,60  ft). 
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Much  of  the  overtopping  data  obtained  during  these  tests  has  been 
presented  by  Saville  (2)  and  by  Saville  and  Caldwell  (l)  and  summarized 
in  CERC's  TR-U  (3);  however,  in  this  latter  publication  the  data  were 
presented  in  dimensional  form  making  their  general  application  difficult. 
In  keeping  with  the  decision  to  present  information  in  CERC's  Shore 
Protection  Manual  (ii)  in  dimensionless  form  whenever  practicable,  the 
overtopping  data  was  reanalysed  and  an  empirical  expression  derived. 

The  broad  range  of  model  scales  used  in  the  overtopping  experiments 
also  provide  an  opportunity  to  investigate  the  effect  of  model  scale 
on  test  results.  A summeiry  of  overtopping  test  Conditions  investigated 
is  given  on  Table  1. 


Figure  1.  Definition  of  terms. 


« 


Numbers  in  parentheses  correspond  to  references  listed  in  Appendix  I. 
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VIBLE  1 - SUMMARY  OP  OVERTOPPING  TEST  CONDITIONS 


R«ngt  of 
W«v«  H«l«htt 


tnvittqaf  d 


Smooth  vertical  wall 
1 on  1-1/2 
amooth  alopa 
1 on  3 amooth  alopa 
1 on  1-1/2 
atappad  alopa 
1 on  1-1/2 

riprap  faced  alopa 
Curved  wall 
Recurved  wall 


Plunger 

type 


TABU  2 


AGREEMENT  BETWEEN  MEASURED  AND  CAICUUTED  OVERTOPPING  RATES 
(Ualng  8PM  publlahed  valuaa  of  d.  and  0^.  baaed  on  1 to  17  acale  data) 


Correlation 


of  Polnta 


Riprap  Face 

1 on  1-1/2 


StaK>ad  Pace 
1 an  1-1/2 


on 

3 amooth  alopa 

193 .5  m long 

Bulkhead 

48.8  to  140.2  cm  0.386  10  10.12  eac 

on 

6 amooth  alopa 

4.S7  m wide 
6.10  m deep 

type 
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DIMENSIONAL  AKALYSIS 


The  variables  describing  the  overtopping  of  a given  structiire  are 
depicted  on  figure  1.  They  include: 


* deepwater  wave  height  [L] 

T » wave  period  [T] 

g 3 gravitational  acceleration  [L]  [T] 

2 -1 

Q » overtopping  rate  (volume  per  unit  time  (L)  [T] 

per  unit  crest  length) 

R » runup  height  measured  vertically  from  (L) 


the  still  water  level  (SWL)  (e.g.  the 
height  to  which  the  water  would  runup 
if  the  structure  were  high  enough  to 
preclude  overtopping) 

d^  ■ water  depth  at  the  structure  toe  [LJ 

h = height  of  the  structure  crest  above  [L] 

the  bottom 

2 

V * kinematic  viscosity  (L)  [T] 


8 ® structure  slope 


(dimensionless] 


plus  any  other  geometric  parameters  necessary  to  describe  the  various 
structure  types.  A dimensional  analysis  of  the  preceding  9 variables 
having  2 dimensions  gives  the  following  dimensionless  terms: 


* relative  water  depth  at  the  structure  toe 
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F * 


Q * 


R * 


e 


h’/6t2 

r/h' 

o 

Q"/gH;3 

0 

' 2 

H 

o 

vT 


= wave  steepness  parameter 

® relative  height  of  structure  crest  above  SWL 

= relative  runup  or  height  of  structure  crest  required 
to  preclude  overtopping 

= relative  overtopping  rate 

= structure  slope,  and 

* a Reynolds'  number. 


The  phenomenon  is  scaled  primarily  according  to  Froude  similarity; 

however,  the  Reynolds*  number  serves  as  a measure  of  any  scale  effects. 

Other  formulations  of  R are  possible,  the  present  one  having  been 
e 

adopted  for  i^s  simplicity. 

Generally  it  is  not  permissible  to  eliminate  dimensionless  terms  by 

combining  them  unless  an  analytic  or  empiriccil  relationship  between  two 

of  the  variables  is  known.  If  it  is  assumed  that  such  a satisfactory 

relationship  is  available  for  the  runup  R,  the  overtopping  rate  can  be 

expressed  in  terms  of  R and  the  ratio  F/F  = (h-d  )/R  can  be  substituted 

o s 

for  F and  F^.  The  preceding  dimensionless  terms  are  obvioxisly  not  the 
only  combinations  of  terms  possible;  however,  they  were  selected  after 
considerable  trial  and  error  because  they  provided  the  greatest  possibility 

for  keeping  dimensionless  variables  constant  and  investigating  the  varia- 

« 

tion  of  Q with  individu€d.  parameters, 

DATA  ANALYSIS 


For  a given  structure  and  set  of  incident  wave  conditions  (e.g. 

constant  H^/gT  and  0),  the  dimensionless  overtopping  rate,  Q 

was  plotted  against  the  dimensionless  crest  height,  F/F  * (h-d  )/R. 

o s 


► 


L 
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A typical  plot  shoving  tvo  data  sets  differing  only  in  model  scale,  is 
shown  in  figure  2.  Generally,  all  data  sets  ^en  plotted  semi-logarith^ 


mically  exhibited  a linear  variation  of  Q with  F/F  for  small  values  of 

• o 


F/F^;  €il80,  the  value  of  Q must  approach  zero  as  the  relative  crest 


height,  F/F^  approaches  1,0  (i.e.,  as  the  crest  of  the  structure  approaches 


the  limit  of  wave  runup).  The  curve  therefore  approaches  F/F  * 1.0  asym- 

o 


ptotically  on  the  semi-logarithmic  plot.  The  hyperbolic  tangent  function 
exhibits  identical  behavior;  hence,  an  equation  of  the  form. 


jT  * a tanh 
o 


log 


(1) 


was  used  to  approximate  the  data.  Here  a and  are  empirical  coefficients 


to  be  established  by  comparing  the  equation  with  the  data.  The  value  of 

a generally  establishes  the  shape  of  the  curve  since  it  is  the  slope  of  the 

• * 

curve  at  F/F  * 0.  Q represents  the  value  of  Q for  a structure  with  its 
o o 

crest  elevation  at  the  SWL.  Figure  3 depicts  equation  1 for  various  values 
of  o.  To  establish  values  of  a and  Q a transparent  template  was  made  of 


figure  3 and  used  as  an  overlay  to  Q vs  E/F^  data  plotted  at  the  same  scale 


on  semi-logarithmic  graph  paper.  By  moving  the  template  vertically  until 
one  of  the  curves  coincided  with  the  trend  of  the  data,  the  value  of  a could 


be  directly  determined.  The  value  of  Q was  determined  by  reading  the  value 

• o 


of  Q where  the  o curves  intersected  F/F  = h-d  /R 

o s 


0.0  on  the  data  plot. 


Thus,  by  overlaying  the  template  to  each  data  set,  values  of  a and  Q were 

o 


established  for  each  structure  type  and  set  of  incident  wave  conditions. 
Interestingly,  the  form  of  equation  1 is  such  that  it  could  be  used  to 
describe  the  overtopping  of  all  of  the  structures  for  which  data  were 
available;  consequently,  figures  similar  to  figure  **  could  be  prepared 


for  each  structure  type.  Such  figures,  which  give  a and  Q as  functions 


of  ^ given  structure  type,  are  presented  in  the 


SPM  (U)  for  other  structure  slopes  and  types. 


data  analysis,  1 on  3 smooth  slope. 
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3y  substituting  the  dimensionless  variables  into  equation  1 and 
solving  for  Q,  one  finds. 


equivalently,  since  teuih  ^ (^)  = ^ l°g^  • 


R+h-d 

( 

'R-h+d  ’ 


(3) 


Either  equation  2 or  3 can  be  used  in  conjunction  with  figures  such  as 
figure  U to  determine  overtopping  rates. 

To  evaluate  the  ability  of  equation  2 or  3 to  predict  the  overtopping 

« 

rates  measured  in  the  experiments,  the  values  of  a and  as  published  in 
the  SPM,  were  used  with  equation  2 and  computed  overtopping  values  compared 
with  measured  values.  Table  2 presents  the  correlation  coefficients  found  in 
the  analysis.  In  general,  agreement  was  excellent;  the  worst  case  was  for 
the  vertical  wall  data  with  r * O.98O.  (The  small  number  of  data  points 
for  the  recurved  wall  make  the  correlation  analysis  for  that  structure 
inconclusive). 

Subsequent  to  publication  of  the  SPM,  further  analysis  of  the  data 

« 

for  smooth  slopes  was  undertaken  in  an  attempt  to  relate  a and  to 

incident  wave  conditions  and  structure  slope.  For  a given  slope,  the 

variability  of  a with  incident  wave  conditions  was  relatively  small, 

suggesting  that  an  average  a could  be  used  and  the  data  reanalysed 

to  establish  the  Q value  that  best  fit  the  data  for  the  average  o. 
o 

The  average  value,  a,  is  shown  on  figure  5 for  four  smooth  structure 
slopes  with  data  obtained  at  three  different  scales. 

The  effect  of  decreasing  model  sccde  seems  to  result  in  a more  rapid 
drop-off  of  the  overtopping  rate  with  increasing  structure  height,  (see 
figure  2).  This  effect  is  also  related  to  the  value  of  R used  to  compute 
F/F^,  however,  an  expression  relating  a with  structure  slope  (smooth 
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Equation  7 is  plotted  on  figure  T for  both  t » l/(2w)  (labelled  linear 
theory)  and  for  e values  determined  from  cnoidal  wave  theory.  Also 
shown  on  figure  7 are  the  values  determined  by  the  analysis  for 
constant  a described  above.  In  general,  data  from  the  1 to  17  scale 
tests  falls  well  below  the  cnoidal  wave  curves  for  all  slopes  for  which 
data  were  available.  Data  from  the  1 to  2 1/2  scale  tests  (squares  on 
figure  7),  however,  are  in  general  conformance  with  the  cnoidal  curves 

I 

when  d /H  >1.5  (non  bre€dcing  waves).  Agreement  appears  best  for  long 
s o 

waves  of  small  steepness.  Steeper  waves,  and  waves  that  break  have 

# 

considerably  lower  values  of  than  predicted  by  the  linear  theory 
or  cnoidal  theory  curves. 


SCALE  EFFECTS 


Scale  effects  arise  because  of  an  inability  to  model  all  aspects 
of  a phenomenon  simultaneously,  usually  because  of  the  limited  range  of 
fluid  properties  practically  achievable;  hence,  to  achieve  both  Froude 
and  Reynolds'  similarity  in  a model,  the  model  fluid  viscosity  would  have 
to  vary  as  the  scale  ratio  to  the  3/2  power  times  the  prototype  fluid 
viscosity.  Therefore,  if  as  is  usually  the  case,  the  same  fluid  is  used 
in  both  model  and  prototype,  only  one  similirity  lav  can  be  satisfied  at 
the  expense  of  the  other.  If  surface  tension  is  also  a factor,  the  problem 
is  even  more  complex.  Since  the  wave  overtopping  phenomenon  is  dominated 
by  wave  motion,  a Froude  modelling  law  governs;  however,  turbulent  and 
viscous  dissipation  also  influence  overtopping;  thus  Reynolds'  effects 
must  also  influence  the  phenomenon.  If  viscosity  is  included  in  the 
original  dimensional  analysis,  a Reynolds'  number  arises  from  the  analysis 
and  serves  as  a "scale  factor"  or  as  a measure  of  scale  effects.  The  for- 
mulation of  the  Reynolds*  number  adopted  to  investigate  scale  effects  was, 

R * H ^/\<r. 
e o 

The  influence  of  R on  an  individual  data  set  for  which  all  other 
e 

dimensionless  variables  are  constant  is  shown  in  figure  2.  Both  sets  of 
points  on  the  figure  are  from  what  were  termed  1 to  IT  scale  tests;  how- 
ever, the  wave  height  in  the  one  test  was  10.76  cm  (solid  circles,  R * 

e 
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8030)  while  in  the  other  test,  the  wave  height  was  21,52  cm  (open  circles, 

F = 22700),  The  greatest  influence  of  scale  appears  when  the  structure 
e 

crest  is  near  the  limit  of  wave  runup  (e.g.  at  the  tail  of  the  curve  near 

F/F  = (h-d  )/R  •*  1,0),  It  is  in  this  relatively  thin  layer  of  fluid  that 
o s 

viscous  effects  can  be  expected  to  increase  in  importance.  Inasmuch  as 
it  is  the  shape  of  the  curve  that  is  affected,  it  is  the  value  of  a that  is 
to  a limited  extent,  influenced  by  scale.  This  slight  dependence  of  a 

on  scale  is  shown  on  figure  8,  For  the  1 on  3 slope,  the  SPM  values  of 

* ' 2 

a are  plotted  against  R for  a broad  range  of  d /H  and  H /gT  values . 

e s 0 o 

While  it  is  difficult  to  draw  firm  conclusions  regarding  scale  effects 
I 1 2 

from  figure  8 since  d /H  and  H /gT  varied,  there  appears  to  be  a decrease 
so  o 

in  a with  increasing  R . Note  that  the  test  designated  as  the  1 to  30 
6 1 2 

scale  (triangular  points)  when  characterized  by  /vT  as  a ’’scale 

factor”  are  actually  not  at  a smaller  scale  than  the  1 to  17  scale  tests 
(circles)  since  the  smaller  wave  heights  also  had  correspondingly  smaller 
wave  periods, 

« 

The  effect  of  scale  on  Q appears  to  be  negligible  since  no  systematic 

# ® # 

variation  of  Q with  R could  be  found.  Since  Q represents  the  dimension- 
o e o 

less  overtopping  rate  for  a structure  with  its  crest  at  the  SVL,  a relative 

ly  thick  layer  of  flow,  insensitive  to  viscous  effects,  tends  to  minimize 

« 

any  dependence  of  on  for  smooth  slopes. 

In  general,  it  appears  that  the  primary  influence  of  scale  on  the 
overtopping  of  smooth  sloped  structures  is  through  scale  effects  manifest 
in  the  runup  phenomenon.  In  the  data  analysis  described  herein,  the  actual 
measured  or  "correct"  runup  was  used;  consequently,  the  influence  of 
scale  effects  on  the  runup  was  not  considered.  Other  investigators  have 
demonstrated  that  small  scale  runup  experiments  tend  to  significantly 
underpredict  runup  at  Isirger  scales.  The  observed  effect  of  scale  on 
overtopping  itself  appears  to  be  relatively  small;  thus  if  the  runup 

values  used  in  equations  2 or  3 are  corrected  for  scale  effects,  the  pre- 

« 

dieted  overtopping  rates  (using  appropriate  a and  values)  will  provide 
good  estimates  of  prototype  overtopping  rates. 
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A levee  is  constructed  with  a 1 on  3 slope  to  Just  preclude  over- 
topping by  a wave  1.05  m high  when  the  water  depth  at  the  structure  toe 
is  1.5  m deep.  If  the  design  wave  period  is  5 sec,  determine  the  rate 
at  which  water  is  carried  over  the  structure  by  the  design  wave  if  a 
storm  surge  of  0.6  m occurs.  The  unrefracted  deepwater  wave  height  can 

be  estimated  from  the  linear  theory;  H/H  = 1.073  for  d /L  = 0.038^*; 

I o s o 

hence,  H = 1.05/1.073  = 0.98  m.  From  reference  the  relative  runup 

* ° t 

is  R/H  = 2.2,  which,  when  corrected  for  scale  effects,  becomes  R/H  = 
o o 

2.2(1.12)  = 2.U6U,  Thus,  R = 2.^6U(0.98)  = 2.Ul  m and  the  height  of 
the  structure  crest  above  the  bottom  is  h = I.5  + 2.Ul  = 3.9I  m.  To 

I 

determine  the  overtopping  rate,  calculate  d /H  = (I.5  + .6)/0.98  = 

10  p so 

2.IU  €uid  H /gT  = 0.98/9.81(5)  = O.OOUO.  From  figxire  I4 , interpolating, 

° « 

o * 0.08  and  = 0.02U.  From  equation  2, 

a = [(9.8l)(0.02lt)(.98)^j  exp  tanh"^  } 

Q = O.U71  exp  |-2.713  tanh"^  (0.75l)j 
Q = 0.U7I  exp  |-2.713  (.975)J 
Q = 0.033^*  m^/sec-ra 

If  the  average  a value  (equation  1*)  had  been  used  with  the  cnoidal  curves 
of  figure  7,  the  value  of  a would  have  been  0.0765  and  S*  = 0.033.  Using 

3 o 

; these  values  in  equation  2 gives  Q * 0.03l»7  m /sec-m. 

CONCLUSIONS 

' 1 ^ t r - /«  - Kt  f * . 

1.  Equation  2 or  3 along  with  the  empirically  established  values  of  a and 
of  reference  .^can  be  used  to  predict  overtopping  rates  for  various 
structure  slopes  and  structure  types  if  accurate  predictions  of  runup  are 
available. 
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2.  Alternatively,  for  smooth  slopes,  if  conservative  (high)  estimates 

are  desired,  eq.uation  2 or  3 can  he  used  with  the  average  a given  hy 

« 

equation  and  the  value  of  Q given  hy  equation  7 with  e determined 
from  figure  6.  This  corresponds  to  using  the  values  given  hy  the 
cnoidal  cxirves  of  figure  7. 

3.  Scale  effects  in  overtopping  tests  arise  primarily  in  modifying 
the  runup.  If  scale  effect  corrected  values  of  runup  are  used  in 
equations  2 or  3,  the  predicted  overtopping  values  will  have  been 
corrected  for  most  viscosity  induced  scale  effects. 
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was  obtained  where 

Q = overtopping  rate  (volume/unit  time) 

g = gravitational  acceleration 

= unrefracted,  deepwater  wave  height 

h = height  of  the  structure  crest  above  the  bottom 

d^  = water  depth  at  the  structure 

R = the  vertical  height  of  runup  above  the  SWL  (if  the  structure  were 
high  enough  to  preclude  overtopping),  and 

Q*,  a = empirical  coefficients  that  depend  on  incident  wave  steepness, 
and  the  water  depth  to  wave  height  ratio  d^/H^. 

Q*  is  a dimensionless  measure  of  the  overtopping  rate  when  a structure 
crest°is  at  the  SWL  elevation  (e.g.,  zero  freeboard  or  h - d^  = 0) . The 
dimensionless  coefficient,  a,  is  a measure  of  the  rate  at  which  overtopping 
decreases  with  increasing  freeboard.  It  varies  primarily  with  structure 
slope  and  shape. 

Surprisingly,  the  equation'  with  appropriate  values  of  a and  Q*  was 
found  to  be  valid  for  all  structure  geometries  for  which  data  were  analyzed. 
Data  were  available  for  smooth,  impermeable  slopes  (from  vertical  to  a 1 on 
6 slope),  stepped  and  rubble-covered  slopes,  and  curved  and  recurved  sea- 
walls. Comparisons  between  predicted  and  measured  overtopping  rates  show 
excellent  correlation  with  the  worst  case  (vertical  wall)  correlation 
coefficient  being  r = 0.980. 
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